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Abstract 



We are using computer models to explore the observational sensitivity to changes in atmo- 
spheric and surface properties, and the detectability of bio signatures, in the globally averaged 
spectra and light- curves of the Earth. Using AIRS (Atmospheric Infrared Sounder) data, 
as input for atmospheric and surface properties, we have generated spatially resolved high- 
resolution synthetic spectra using the SMART radiative transfer model, for a variety of con- 
ditions, from the UV to the far-IR (beyond the range of current Earth-based satellite data). 
We have then averaged over the visible disk for a number of different viewing geometries to 
quantify the sensitivity to surface types and atmospheric features as a function of viewing 
geometry, and spatial and spectral resolution. 

According to our model, in the case of Earth several atmospheric species can be identified 
in disk- averaged spectra and potentially detected depending on the wavelength range and re- 
solving power of the instrument. O3, H2O, O2, CO2 (and probably oxygen dimer (02)2) are 
clearly visible in the optical part of the spectrum, CO2, O3, H2O in the MID-IR (CH4 and 
N2O are hard to discriminate because obscured by the strong water absorption). 
The detectability of the red edge (a distinctive signature of vegetation) is affected by: 1) the 
extent of the vegetation over the planetary surface, 2) how much the vegetation is obscured 
by clouds, and 3) the strength of the signal for the particular vegetation type. According to 
our model, the land surface cover of vegetation on Earth seems to be adequate to produce a 
disk- averaged signal that is strong enough to show the red edge in disk- averaged spectra, even 
when the signal is averaged over the daily time scale. To quantify how much the vegetation is 
obscured by clouds, we selected some specific bands as indicators of photo synthetic activity. 
The possibility of detecting phyto -plankton appears to be more challenging. 
These results have been processed with an instrument simulator to improve our understand- 
ing of the detectable characteristics of Earth-like planets as viewed by the first genera- 
tion extrasolar terrestrial planet detection and characterization missions ( Terrestrial Planet 
Finder/Darwin and Life finder). The wavelength range of our results are modelled over 
are applicable to both the proposed visible coronograph and mid-infrared interferometer TPF 
architectures. We have validated this model against disk- averaged observations by the Mars 
Global Surveyor Thermal Emission Spectrometer (MGS TES). This model was also used to 
analyze Earth-shine data for detectability of planetary characteristics and biosignatures in 
disk- averaged spectra. 



1 Introduction 



The principal goal of the NASA Terrestrial Planet Finder (TPF) and ESA Darwin 
mission concepts is to detect and characterize extrasolar terrestrial planets. TPF-C 
(coronograph) , TPF-I (interferometer) and Darwin are expected to survey nearby 
stars and directly detect planetary systems that include terrestrial-sized planets in 
their habitable zones (TPF website, TPF book, Darwin website). 

This first generation of terrestrial planet exploration missions will only be able to 
resolve the planets as point sources. The distance and the relatively low brightness of 
the terrestrial sized planet compared to its parent star pose a daunting technological 
challenge: the spectral information provided will be necessarily an average over the 
visible disk and also over the exposure time, which will probably be non negligible 
(hours or days depending on the target). 

The task for the scientific community will be to interpret this space-time-averaged 
single pixel, in which the information of the whole planet is irreversibly collapsed. 
Due to the nature of the problem, in the most general case the interpretation of the 
recorded spectrum will not be unique. A family of solutions will provide an equally 
good explanation of the spectral features observed. This expected degeneracy is due 
in part to a lack of spectral, spatial and temporal sensitivity of these averaged spectra 
simulating the observations. Another source of degeneracy is given by space-time 
symmetries: those can be studied and predicted with our model as well. 

Previous research on detectability of planetary characteristics in the disk-average, 
has focused on measurements of earthshine on the unilluminated side of the moon 
(Woolf et al., 2002; Arnold et al., 2002; Goode, Montahes Rodriguez, Palle et al., 
2001-2003), and models of diurnal photometric variability on an Earth-like planet 
(Ford et al., 2001). 

Ford et al. (2002) developed a model performing Monte Carlo integrations, with 
single scattering over a planet with different surface types specified by a map. They 
used this code to evaluate the photometric variability in some specific bands in 
the optical due to daily rotation and seasons. Although providing a good estimate 
of surface reflectivity in the optical, their model did not include an atmosphere or 
realistic clouds: both phenomena can significantly change the reflectivity of a planet. 

Goode, Montahes Rodriguez, Palle et al., (2001, 2003) observed for several years 
earthshine reflected by the moon to determine precisely the global and absolutely 
calibrated bond albedo of the Earth, and to study atmospheric phenomena govern- 
ing the Earth's radiation budget. Woolf et al. (2002) examined earthshine data 
and speculated on the roles of different atmospheric and surface features (clouds, 
ozone, water, oxygen, Rayleigh scattering, detectability of the vegetation signal). 
Arnold et al. (2002) used a collection of earthshine data (400-800 nm) to derive nor- 
malized Earth albedo spectra and extract the signature of Earth vegetation. These 
observations, although extremely valuable, are necessary limited to views of the 
Earth available from astronomical observatories and over a relatively short time 
span. There are viewing geometries that can not be obtained via earthshine (eg. 
South or North Pole views) and that might be of interest for missions like TPF 
and Darwin, which may observe planetary systems at significant inclinations. The 
advantage of a comprehensive theoretical model that realistically includes the contri- 
butions of surface atmosphere and clouds is the possibility to generate the full range 
of scenarios (clouds, atmospheric gases, geometries) that are not available from a 
single observational data set. 

To better understand these issues, we built a model which reproduces disk- 
averaged spectra as faithfully as possible, starting from a description of a spatially 
resolved terrestrial planet on diurnal or seasonal time scales. The model we imple- 
mented, is a useful tool for exploring the sensitivity of the time-disk-averaged spectra 
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to local properties and processes on time scales shorter than the integration time of 
the instrument. These results were then processed with an instrument simulator to 
improve our understanding of the detectable characteristics as viewed by the first 
generation of extrasolar terrestrial planet detection and characterization missions. 

We present here the first simulations of disk-averaged spectra and relative light 
curves (obtained by integrating the spectral signal over a specific band) of Earth 
from the Infrared to the visible. This paper is the continuation of a previous one 
about Mars (a good example of a likely abiotic planet). We refer to that paper for 
the details of the model and for a useful planetary comparison (Tinetti et al., 2004). 

Disk- aver aged spectra show sensitiveness to seasons, viewing geometries, phases 
and surface signatures. The red edge, in particular, is a distinctive signature of 
plants. Astrobiologists are interested whether this biosignature could arise on an 
extrasolar planet and be detectable with TPF. The question of the origin of the red 
edge is beyond the scope of this paper (see discussion section), but, assuming that 
it is plausible, our model and the earthshine data can help to assess the challenges 
for detecting it in a disk-averaged spectrum. The primary challenges are: the spa- 
tial extent of vegetation over the planetary surface, and how much they might be 
obscured by clouds. Large part of the result and discussion sections is dedicated to 
the selection of the bands to be used as indicators of photosynthetic activity. Some 
of them lose significantly their sensitivity when clouds are included in the simula- 
tions. The choice of the bands needs necessarily to take into account those effects. 
Clouds in fact deeply affect the intensities and shape of the Earth spectra, masking 
or reducing in some cases or enhancing (depending on cloud altitude, composition 
and optical depth) in some other cases some atmospheric and surface features. This 
consideration can probably be extended to the case of a general terrestrial planet, 
even if the specific results valid for the case of Earth and water clouds might not be 
applicable to the case of exotic condensate clouds and atmospheres. 

The model has been validated against disk-averaged observations of the Earth 
by the Mars Global Surveyor Thermal Emission Spectrometer (MGS TES, fig. EJ) 
(Christensen, Pearl, 1997) and earthshine spectra (fig. IH|) (Woolf et al. 2002). The 
model generates a variety of products including disk averaged synthetic spectra, 
light-curves and the spectral variability at visible and IR wavelengths as a function 
of the viewing angle. These tools have been used for simulations of an increasingly 
cloudy /forested/oceanic Earth (fig. IH1 and [T5|) and analyzed to determine the de- 
tectability of biosignatures on an Earth-like planet (e.g. red-edge signal, fig. El El 
and El presence of plankton, fig. fTHIand fTTj) . We have therefore examined the sensi- 
tivity of the red-edge detect ability to varying 1) planetary geometric views (varying 
visible land cover), 2) clouds (effect of both cloud type and amount of cover) over the 
current Earth's surface, and 3) the concentration of chlorophyll in phytoplankton. 
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2 Modeling approach 



The core of the model is a spectrum-resolving (line- by-line) atmospheric/surface 
radiative transfer model, SMART, (Meadows and Crisp, 1996; Crisp, 1997). AIRS 
(Atmospheric Infrared Sounder) data were used as input for Earth surface albedos 
and atmospheric properties (Fishbein et al., 2003). 

These data were processed and mapped on the sphere using a pixelization scheme 
created for Planck and WMAP missions (Healpix) (Gorski et al., 1998) (fig. GJ. For 
each of the pixels SMART was run to generate a database of synthetic spectra for 
a variety of viewing angles, illuminations, surface types and cloud cover. Finally, 
a model was created in which the user specifies spatial resolution and observing 
position, and the program selects the appropriate synthetic spectra and creates the 
final view, which can then be averaged on the disk (Tinetti et al, 2004). 

In the following paragraphs we give a fairly detailed description of the models, of 
the data we used as input to model, and of the methods we have selected to estimate 
some specific parameters or indexes. 



Visible 




Figure 1: For producing disk- averages that could be viewed from any direction, spherical coordi- 
nates were not the best choice: parallels and meridians do not divide the sphere into equal area tiles, 
the pole is a singularity, the geometry is not invariant for a general rotation. 

We chose to use Healpix (Hierarchical Equal Area and iso-latitude Pixelization) grid as tiling pro- 
gram for our model. Healpix, originally designed for the ESA Planck and NASA WMAP missions 
(Gorski et al., 1998), is a curvilinear partition of the sphere into exactly equal area quadrilaterals 
. This peculiarity is particularly desirable since the geometry does not considerably change when we 
look at the planet from different viewing angles, and the contribution of each region of the sphere is 
accounted for with the same weight. The ability to easily degrade and upgrade the resolution of the 
map, is also useful. 

We resolved the Earth with 48 pixels for mapping the atmospheric properties and more than three 
thousands pixels for tiling the earth's surface. 
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2.1 Input data 



Atmospheric properties The atmospheric temperature, water vapor trace-gas 
and cloud distributions needed to construct the model are extracted from the NCEP 
Global Forecast System (NCEP, 1988; Kanamitsu et al, 1991; Kanamitsu, 1989) and 
supplemented by the Upper Atmospheric Research Satellite (UARS) and Harvard 
University tropospheric ozone climatologies by the AIRS Level 2 Simulation System 
(Fishbein et al., 2003). The NCEP forecasts are produced every three hours and have 
l°times l°horizontal resolution and extend upto lOhPa. The UARS climatologies 
are monthly zonal means with 10°latitude resolution and extent above 0.01 hPa while 
the Harvard tropospheric ozone climatology (Logan, 1998) is an annual mean with 
4°times 5°resolution. 

The simulation system samples the fields as the AIRS instrument would view 
the earth from its 705 km altitude sun-synchronous 1:30 PM local-time ascending 
equator-crossing polar orbit. The satellite orbits the earth 14.6 times per day col- 
lecting 3 million observations daily. The AIRS scan produces gores between the 
orbit passes near the equator (Fig. which have been filled using a Gaussian inter- 
polation. The data have been remapped to a spherical grid, 37 by 72 (latitude by 
longitude) (Fig. EJ). 




Figure 2: Left: This figure shows the actual AIRS solar albedo map. The dark blue strips are 
the regions where there are no data points. Right: This figure shows the same data after using the 
two-dimensional Gaussian interpolation. The strips are no longer visible but the contours are less 
defined. 

The spatial interpolation method assigned each of the 3 million data points (taken along the near 
polar orbit path) to a regular latitude-longitude grid by performing a weighted average, using a two- 
dimensional Gaussian distribution centered at the location of the data point. The average value 
at each grid point was normalized by summing the contributions of 3 million Gaussians at each 
latitude/longitude grid point. This procedure was repeated for each parameter - temperature and 
trace gases. 

Solar spectrum The upper boundary condition of the model is specified by the 
downward solar flux at the top of the atmosphere. The solar spectrum used is a 
high-spectral-resolution spectrum compiled from satellite observations at UV wave- 
lengths (Atlas-2 SUSIM) and stellar atmosphere models at longer wavelengths (Ku- 
rucz, 1995). The position of the sun in our simulation was derived from the Jet 
Propulsion Laboratory Horizons Ephemeris System (Giorgini et al., 1997). 



Molecular Line lists The optical properties of the atmosphere were simulated 
using molecular line databases to generate the line-by-line description of infrared 
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Fi glire 3: Examples of temperature-pressure (fig. on the left) and gas-mixing ratios profiles 
(H2O, O2, O3, CH4, CO2, CO, N2O, fig. on the right) used for our model 

vibration-rotation bands of gases, and UV and visible cross-sections of gases associ- 
ated with electronic and pre-dissociation transitions (e.g. ozone). 

For our model, we have used the HITRAN2K (High resolution Transmission 2K) 
database [Rothman et al., 2003]. 

Albedo/Surface Types The land cover map provided by AIRS assigns to each 
point one of eighteen possible surface types (Fig. 0J. The contribution of each 
material is determined by land fraction, the amount of vegetation, and the types 
of vegetation defined by the International Geosphere Biosphere Programme (IGBP) 
land use surface classification (Loveland and Belward, 1997). For our purposes, 
these 18 divisions were then reduced to six general surface types: ocean, forest, 
grass, tundra, ground, and ice. 

These different surface types were simulated using as input to the radiative 
transfer model six measured, wavelength dependent, surface reflectance spectra, 
a k (X), k = 1...6. 

The IGBP map (Fig. 0J) and the solar albedo map (fig. Ej) must be consistent. 
The surface reflectance spectra a^(A) must then be recalibrated in order to give 
finally the same global average albedo provided by the AIRS albedo map. The 
scaling factor for each surface type can be calculated by comparing the solar albedo 
value from the AIRS map (fig. 0|) with the numerical integral in the 0.49-0.94 /im 
range (the wavelength range of the AIRS instrument which provided the albedo map) 
of the a^(A). Using the simplified land cover map (6 surface types only), for each 
of the 48 atmospheric pixels, the fractional area covered by each of the six surface 
types was determined. The radiative transfer code was run only for the surface types 
with a non-negligible fractional area (greater than 5% of the entire pixel area). 

Clouds Clouds vary considerably in amount and type over the globe. They have 
very important effects on longwave and solar energy transfer in the atmosphere, as 
we can see from Table 1 (Harrison et al. 1990), (D. Hartmann, 1994), (Oort and 
Peixoto, 1992). 

For radiative calculations optical and radiative properties must be specified. Nor- 
mally, water clouds have relatively weak solar absorption, but they effectively scatter 
solar radiation back toward space and so have high reflectivities. Thick clouds be- 
have almost like black bodies for long wave radiation: they absorb all incident long 
wave radiation and emit like a blackbody with the temperature of the atmosphere 
at that the same level of the cloud. A simple approach is to specify the properties 
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Fi gure 5: Reflectance in the optical/ NIR for the six surface types we used in our model 



of three types of clouds: high (cirrus), medium (cumulus) and low (stratus). 



Type SW reflectivity SW absorptivity % area 

High (cirrus) 0.21 ' 0.005 0.228 

Medium (cumulus) 0.48 0.02 0.09 

Low (stratus) 0.69 0.035 0.313 
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Table 1 - Values of cloud shortwave reflectivity and absorptivity and fractional area coverage 
assumed in (Manabe and Strickler, 1964) 

In our model, the wavelength dependent optical properties of liquid water droplets 
were computed using Mie scattering theory (Crisp, 1997), (Liou, 2002). Cirrus clouds 
were parameterized as polydispersions of hexagonal crystals, and their optical prop- 
erties were derived using geometric optics (Muinonen et al., 1989). The liquid water 
and ice refractive indexes were obtained from (Segelstein,1981) and (Warren, 1984) 
respectively. 

The top and bottom pressures and the optical depths were taken from AIRS 
Level 2 Simulation System (Fishbein et al., 2003). Those sets of data are sometimes 
underestimating the presence of low clouds. We have used both AIRS and MODIS 
data (MODIS website) for the cloud maps. 
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2.2 Models 



Most of the models listed here were already described in a previous paper about 
Mars (Tinetti et al., 2004). We refer to that work for a more detailed explanation 
and list of references. 

LBLABC LBLABC is a line-by-line model (Meadows and Crisp, 1996) that gen- 
erates monochromatic gas absorption coefficients from molecular line lists, for each 
of the gases present in the atmosphere at a specified temperature and pressure. This 
model was designed to evaluate absorption coefficients at all atmospheric levels over 
a very broad range of pressures (10 -5 to 100 bars), temperatures (130 to 750 K), 
and line center distances (10 -3 to 10 3 cm -1 ). 

SMART SMART is a multi-stream, multi-level, spectrum-resolving (line-by-line) 
multiple-scattering algorithm designed to generate the high-resolution synthetic spec- 
tra of planetary atmospheres (Meadows and Crisp, 1996). Using a high-resolution 
spectral grid, it completely resolves the wavelength dependence of all atmospheric 
constituents, including absorbing gases (infrared absorption bands, UV predissoci- 
ation bands, and electronic bands) and airborne particles (clouds, aerosols) at all 
levels of the atmosphere, as well as the wavelength- dependent albedo of the planet's 
surface, and the spectrum of the incident stellar source. 




Figure 6: This figure shows the three input angles for SMART: the solar zenith angle (0i), 
the viewing angle zenith angle (0 r ) and the viewing azimuth angle (<p r — <\>%)- The solar azimuth 
angle (<j)i) has been chosen as the reference direction. SMART was run for each pixel for nine solar 
zenith angles (0, 15, 30, 45, 60, 75, 80, 85, and 90 degrees), four viewing zenith angles (the four 
Gaussian angles: 21.48, 41. 93, 70.73, and 86.01 degrees), and seven viewing azimuth angles (0, 
30, 60, 90, 120, 150, and 180 degrees), six surface types (tundra, ocean, ice, ground, grass and 
forest) and three cloud types (Cirrus, Strato- Cumulus and Alto- Stratus). This library was used for 
the creation of three-dimensional data cubes (two dimensions spatially and one dimension for the 
spectra). In the spatial dimensions, only the pixels that are detectable by the remote sensor will be 
considered, therefore the data cube will be comprised of pixels up to 90 degrees in each direction from 
the sub-sensor position on the planet. The resultant data cube will display one half of the planet's 
surface with the corresponding spectra at each of the detectable pixels. 
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Averaging on the disk We used a very high resolution of the Healpix base (more 
than three thousands) to assign to each of the 48 atmospheric pixels the fraction 
detectable by the sensor, the fration illuminated and the fraction covered by different 
surface types and clouds (we took in consideration 3 cloud types: Strato-Cumulus, 
Alto-Stratus, Cirrus) according to the AIRS-IGBP land cover and cloud maps (AIRS 
-level 2, and MODIS). 

SMART was run for a variety of surface types (six in total), cloud types (3 main 
types), solar zenith angles (nine), viewing zenith angles (four), viewing azimuth 
angles (seven) (fig. EJ). The radiance values for arbitrary angles can be estimated 
using interpolation (linear interpolation for the solar zenith angle and bicubic spline 
interpolation for the viewing angles) 

Using the simplified land cover map and a cloud map, for each of the 48 atmo- 
spheric pixels, the fractional area covered by each of the six surface types and by 
one of three cloud types was determined. We neglected the contribution of surface 
types covering a fractional area of less than 5%. 

The radiance for each pixel can be calculated using the method previously ex- 
plained for the viewing/solar angles and summing the contributions of the all the 
radiances corresponding to different surface and cloud types weighted by their frac- 
tional area. 

The disk-averaged synthetic spectrum can be easily calculated by summing the 
radiance values for all of the detectable pixels at a given wavenumber and dividing 
by the number of pixels. 
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3 Earth: Validation with observed spectra 



We have validated our model against two types of observations. We used MGS-TES 
data taken from Mars Global Surveyor en route to Mars for the IR (Christensen, 
Pearl, 1997) and earthshine data (Woolf et al., 2002) for the optical. 

Earthshine is sunlight reflected by the Earth that can be seen lighting up the dark 
portion of the crescent Moon. To obtain the earthshine, the spectrum of the dark 
Moon is divided by the spectrum of the bright Moon, thereby removing the reflection 
spectrum of the Moon itself and the extra pass through the Earth's atmosphere. The 
result is the disk-averaged outgoing radiation at the top of the Earth's atmosphere. 
This can be directly compared with our model. 

Our model produces an acceptable fit with the IR and optical data as we can see 
in fig. d and El In both these figures we have plotted the recorded (black plot) and 
the synthetic spectra (red plot), and the contribution of the cloud free (yellow plot) 
and totally cloud covered spectra as well (light blue plots). 

In the IR (fig. EJ) the clear-sky spectrum is the uppermost curve (warmest) since 
with no clouds we can observe the warm lower layers in the troposphere and the 
surface. The main effect of adding clouds is to decrease the emitted flux (colder 
temperature where clouds form). The reduced thermal contrast between the cloud 
tops and space also reduces the amplitude of the spectral features: the lowest curve 
(cirrus clouds) is almost a blackbody at the tropopause temperature. The emission 
spike in the middle of the 15 /im CO2 band and 9.6 /im O3 band is due to the Earth's 
stratospheric temperature inversion. H2O absorbs at wavelengths longer than 15 /im 
(rotational bands) and from 5 to 8 /im (vibrational bands). CH4 and N2O are hard 
to detect at (respectively) 6 and 8 /im because their absorption is obscured by the 
strong water absorption. 

In the optical (fig. IH|) the contribution of clouds is more complex: clouds affect 
the atmospheric light paths in several ways compared to the clear-sky conditions. 
Because of multiple scattering inside the clouds, the light path can be either increased 
or decreased, and the absorption of the tropospheric trace gases can be strongly 
enhanced or reduced. Thin high clouds can also reduce the absorption paths below 
the clouds (Wagner et al., 1998). As shown in figure[Hl we confirm the presence of O3 
(Chappuis bands, 0.5-0.7 /im), H2O (~ 0.7, 0.8, 0.9 /im), O2 (B band, 0.69 /im, and 
A band, 0.79 /im) already discussed in (Woolf et al. 2002). Moreover, we believe the 
absorption features ~ 0.6 /im might be explained (C. Miller, private communication) 
by the oxygen dimer (02)2 (absorption bands at 0.577 and 0.63 /im) (Wagner et al., 
2002). The model is unable to reproduce this feature since (02)2 was not initially 
included in our simulations. An extensive discussion about atmospheric species and 
spectral signatures in both IR and visible can be found in (Des Marais et al. 2002). 

The phase and the viewing geometry for the Woolf et al. (2002) optical observa- 
tions were known and were used in our model to reproduce the observing geometry 
precisely. Given the observing geometry and the day of the year, we calculated the 
extreme cases for a completely clear and completely cloudy Earth (fig. [HI yellow 
plot and light blue plots), which must bound the correct fit. However, to model the 
observed spectrum precisely, we used a realistic spatially resolved cloud cover for 
the date of observation (AIRS Level 2 data and MODIS data slightly adjusted to 
obtain a better fit). In doing that, we have tried to reproduce not only the spectral 
shape, but also the absolute brightness, neglecting the possible uncertainties present 
in the measured spectrum (the data were taken through intermittent clouds and 
the method of scattered light subtraction is not precisely known). The best fit was 
obtained considering ^30% of cirrus, ^10% of Strato-Cumulus and ^10% of Alto- 
Stratus clouds spatially distributed. The oxygen A-band at 0.76 /im was particularly 
useful for cloud retrieval (Kuze and Chance, 1994). 



11 







8 10 12 14 16 18 20 

Wavelength (/i.m) 



Figure 7: on £/ie left: disk-averaged spectrum of Earth observed by Mars Thermal Emission 

Spectrometer (TES) (black) (Christensen, Pearl, 1997) compared with the synthetic one produced in 
almost the same conditions with almost 50% of clouds (red) and without clouds (yellow). Plots in 
light blue show the contribution of different cloud-types. Fig. on the right: we show the view of the 
Earth at 17:30 UT on November 24, 1996, from the Mars Global Surveyor (Walker, 1994). 
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Figure 8: Fig. on left: Earth shine spectrum of the Earth observed from the Steward Observatory 
of the University of Arizona (black plot) compared with the synthetic one produced in almost the same 
conditions of illumination and cloud cover (red). The yellow plot at the bottom of the figure, is the 
disk-averaged synthetic cloud-free spectrum. Plots in light blue show the contribution of different 
cloud-types. Fig. on the right: we show the view of the Earth at the time of the measurement (Woolf 
et al, 2002), (Walker, 1994). 
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4 Results 



Effects of Surface types and clouds 

To assess the contributions of the various surface types, we have run some simulations 
covering the Earth's surface varying one environmental variable at a time. The plots 
in fig. [HI show disk- averaged spectra of ocean (dark blue curves), ice (grey curves), 
forest (darker green curves), grass (green curves) and tundra (black curves) worlds. 
The viewing geometry and illumination are the same in all the plots (view: latitude 
0, longitude 0, fully illuminated). In the IR it is impossible to distinguish among 
the different surface types (fig. El top left). This is due to the fact that the surface 
reflectance does not change appreciably in the IR, just the opposite in the optical 
(fig. The light blue plots in the figures show the contribution of clouds - cirrus 
(continuous line), alto-stratus (dashed line) and strato-cumulus (dotted line)-. In 
fig. El top right, we show a zoom into the 0.65-0.85 fim band, where the red-edge 
signal is detectable. Leafy plants reflect sunlight strongly in this band. We will focus 
on the red-edge detectability later on in this section. 

Disk-averaged synthetic spectra of Earth. IR and visible 

Disk-averaged solar and IR spectra of Earth were generated from several vantage 
points (i.e. over the pole, over the equator, etc.) with the same illumination and 
with no clouds (fig. fTU|) . In the optical, the main differences among the spectra are 
due to the changed surface composition. In particular, in fig. E3 we show a day 
in Northern Hemisphere summer where the Southern polar cap is more extended 
(ice has a very high reflectance in the optical, fig. E|) and the Northern hemisphere 
is more forest covered. There are visible discrepancies between the two equatorial 
viewing geometries due to the high presence of land/vegetation in one case and Ocean 
on the other. Both cases simulate a viewing geometry including a large amount of 
vegetation, and the black and light blue plots show a very distinct red edge signal 
starting from 0.7 /im, 

In the IR, on the contrary, the observed variations are due to the horizontal 
temperature gradients. As we expect due to the time of the year, the lowest average 
temperature is the one of the South pole view. 

Sensitivity to phases 

We have run some experiments to explore the sensitivity of disk-averaged spectra 
to different phases. No pronounced differences are seen in the IR, except for those 
caused by changes in the temperature/cloud distribution due to a changed illumina- 
tion (~ 6% for the cloud-free case and ~ 12% for the low-cloud case, compare the 
four plots on the right in fig. fTT|) . Figure [H] shows the results for the optical as well. 
For the cloud free disk-averaged spectra (black plots) there are only slight discrep- 
ancies among the four cases. On the contrary, dramatic changes can be seen for the 
fully cloud covered disk-averaged spectra (light blue cirrus, blue Alto-Stratus, violet 
Strato-cumulus). The differences are both in intensities and shapes. This behavior 
is mostly due to variations in atmospheric pathlengths above the clouds. When the 
visible disk is fully illuminated all the angles from to 90 contribute symmetrically 
to the final integration, when less than only half of the disk is illuminated, only the 
angles close to 90 contribute to the final result. 
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Figure 9: Disk-averaged spectra with single surface type coverage (blue: ocean, grey: ice, yellow: 
desert, green vegetation, black: tundra). Almost no differentiation in the IR (radiance expressed in 
W/m 2 /sr). In the visible (radiance divided by the solar radiation at the top of the atmosphere) on 
the contrary we have pronounced differences in the spectra. The light blue plots in the figures show 
the contribution of clouds: cirrus (continuous line), alto-stratus (dashed line) and strato- cumulus 
(dotted line). In the band 0.7-0.8 \im we can recognize very clearly the red-edge signal (leafy plants 
reflect sunlight strongly in this band, fig. WJf) . 



Light curves 

Time dependent variations in the disk-averaged spectra, or "light curves" can provide 
additional information about spatial variations and periodicities. Fig. [121 left shows 
light-curves for the spectral interval 8-13 /im following the diurnal rotation of the 
planet for three different phases. The phases considered are totally illuminated, 

totally dark and dichotomies. The quantity plotted is — ^ , where Z(A) is the 

disk-averaged radiation, Ai and A2 are the extremes of the chosen interval (8-13 
/im). Although all the three curves clearly show the daily periodicity, the rotational 
variations in this spectral range are of the order of magnitude of 3% for the fully 
illuminated disk and 0.3% for the totally dark disk. It is unlikely that this generation 
of extrasolar terrestrial planets missions will be sensitive to those fluctuations. The 
relatively small horizontal gradients in the thermal are due to the presence of a quite 
thick atmosphere and an Ocean that act as a buffer. On a planet like Mars, where 



14 



0.25 




Wavelength (^m) 




Wavelength (^m) 



Figure 10: Earth solar (fig. at the top, radiance divided by the solar radiation at the top of the 
atmosphere) and IR disk-averaged spectra. In this particular simulation there are no clouds and the 
disks are equally illuminated. We have used AIRS data corresponding to July 20th 2002 as input 
for SMART. Viewing position: black curve, North Pole; light blue curve, South Pole; green curve, 
Equator-longitude 0° , red curve, equator-longitude 180° . 



the atmosphere is one thousand time smaller than the one on Earth and where there 
is not a huge reservoir of liquid water on the surface, heat is exchanged above all 
through conduction and radiation. As a result, the temperature gradients between 
day and night and between different positions on the surface are extreme (Tinetti et 
al, 2004). The fluid component on Earth is playing a fundamental role in the global 
energy balance, by storing heat through the water cycle (subtraction and release of 
latent heat) and in the smoothing of horizontal gradients through convection. 

Fig. [12] right shows light curves in the visible for totally illuminated disks. For 
the chosen bands 0.5-0.6 /im, 0.6-0.7 /im these light curves are almost flat lines. The 
band 0.7-0.8 /im shows on the contrary high variability due to the presence of the 
red edge signal (see next paragraph for a more detailed discussion). 
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Figure 11: The black plots are the cloud free disk-averaged spectra relative to the four different 
phases, the others represent the fully cloud covered disk-averaged spectra (light blue cirrus, blue 
Alto- Stratus, violet Strato- cumulus). 
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Figure 12: Fig. on the left: light-curves following the diurnal rotation of the planet for the 
spectral interval 8-13 fim in the IR for three different phases. Fig. on the right: light-curves following 
the diurnal rotation of the planet for the spectral intervals 0.5-0.6 \im, 0.6-0.7 /jm, and 0.7-0.8 \im 
in the visible. We see almost no variability in the 0.5-0.6 and 0.6-0.7 bands, the 0.7-0.8 band on the 
contrary shows a pattern due to the presence of the red-edge. 



17 



Detectability of surface biosignatures 

Vegetation "Red edge" and NDVI idex The grass and forest albedo curves in 
fig. Elillustrate the red edge, which begins between .680 /im and .700 /im and rises to 
a plateau at .761 /im. The low reflectance in the .450-. 520 /im (blue) region is due to 
absorbance by carotenoids and chlorophylls, and that in the .630-. 700 /im (sometimes 
extending up to .730 /im) region is due only to chlorophyll, with the higher reflectance 
in the green giving the green color of most plants (Tucker, 1978). Because only the 
chlorophylls are actually involved in reduction of molecules relevant to photosynthesis 
(other pigments assist in light harvesting), we are interested primarily in the red 
portion of the visible absorbed by the chlorophylls for detection of photosynthetic 
activity. Some regions of the NIR reflectance plateau may vary according to plant 
physiological status (Pehuelas and Filella, 1998). 

To quantify the strength (hence the potential detectability) of the vegetation 
"red edge", we use a commonly used index of vegetation activity, the Normalized 
Difference Vegetation Index (NDVI), which is a measure of the contrast between the 
red and NIR reflectances. A high NDVI indicates high absorbance by chlorophyll 
and a well-hydrated status. Note that NDVI is only an indicator of photosynthetic 
activity, but its use in any actual quantification of photosynthesis is only correlative, 
because many other vegetation features (e.g. canopy structure, leaf surface charac- 
teristics) can cause variation in NDVI. Nonetheless, it is a useful indicator. NDVI is 
calculated as: 

NIR + R v ; 

where NIR is near-infrared reflectance, and R is red light reflectance. The choice of 
bands for calculating NDVI varies among different workers depending on available 
data. Landsat Thematic Mapper sensors have bands at .63-. 69 /im and .74-1.10 /im 
for the red and NIR, respectively (Tucker, 1978). NOAA uses .58-. 68 /im and .725- 
1.000 /im for the visible and NIR, respectively. MODIS sensors, designed to avoid 
windows of atmospheric absorbance for bands specific to surface remote sensing, 
have the red band at .649-. 679 /im and an NIR band at .855-. 875 /im (Salmonson, 
1990). These ranges avoid atmospheric absorbance by O2 and H2O vapor. 

Here we calculate the disk-averaged NDVI of our model using .65-. 68 /im for the 
red and .75-. 8 /im (band I) or .855-. 875 /im (band II) for the NIR. 

If we calculate the NDVI index for a particular surface type (fig. Ej), we find the 
following results: 
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Table 2 - Values of NDVI indexes (band I and band II) for different surface and cloud types. 
In the clear sky case, only vegetation has a positive NDVI index. 

In a clear sky case the only surface showing a positive NDVI index is vegetation. 
Clouds may produce a positive NDVI index for certain phases, viewing geometries 
and distribution patterns. The numbers we report in table 2 for clouds refer to 
a pessimistic case -totally illuminated, cloud uniformly distributed-. Even if those 
numbers are considerably lower than the vegetation case, clouds may cover a larger 
area of the visible disk compared to vegetation. If we keep 0.03 -band I- and 0.09 
-band II- as thresholds to discriminate unambiguously the presence of vegetation, 
every value between and the thresholds might be explained either by clouds or 
vegetation. 

If we calculate the NDVI index for the case of earthshine discussed in the valida- 
tion paragraph (fig. 8), Woolf et al. data give 0.015 (band I) and 7.-10 -3 (band II). 
Band I seems to indicate that vegetation is present in the disk average, since NDVI 
values for clouds for that particular phase and distribution are all negative. Band II 
is ambigous, since the the NDVI index is positive but below the threshold. 

Red-edge signal and clouds Fig. ITTT1 shows modeled disk-averaged spectra over 
0.6-.95 microns of different views of a cloud-free Earth during a diurnal cycle. The 
NIR reflectance by vegetation is clearly visible, since it is responsible for a variation of 
50% in the signal. The surface cover of vegetation on Earth is evidently adequate to 
produce a strong enough disk-averaged signal to show a red edge feature when there 
are no clouds. Although the initial rise of the red edge at around 680 nm is obscured 
by absorbance by atmospheric water vapor, the latter part of the red edge is captured 
in the rise in the 0.73-0.81 range. The plateau at .745 microns occurs earlier than the 
actual vegetation red edge plateau at 0.76 microns due to atmospheric water vapor. 
The bottom plot in Figure fTTTI shows the sensitivity of NDVI, our red edge strength 
indicator, to different views of the Earth, cloud cover and for two different choices of 
the NIR band: 0.75-0.8 microns (band I) and 0.855-0.875 microns (band II). Here, 
the red edge can not be discriminated when the Pacific Ocean dominates the view, 



19 



while the peak detectability occurs when Africa and the Eurasian continents provide 
a view with the maximum land cover (~ 40%). Considering the values in table 2, we 
can infer that vegetation is clearly present when NDVI index is positive (vegetation 
cover - 20 %). 

Introducing clouds (uniformly distributed) dampens the variability, as we see the 
red edge detectability is progressively reduced for cloud cover of 30% and 60%. For 
the clear sky experiment band II seems to be a better choice than band I, but when 
clouds are taken into account band II produces a more ambiguous signal compared 
to band I (Table 2). What level of NDVI should be accepted as a true signal is a 
subject for the discussion section. 

The surface cover of vegetation on Earth is evidently adequate to produce a 
strong enough disk-averaged signal to show the red edge even when we average over 
the diurnal cycle (fig. fT^j) . This information can be useful for a mission like TPF-C 
which most probably will have to deal with integration times of the order of a day. 

Figure CS1 shows the sensitivity of the disk- averaged and diurnally averaged spec- 
tra to percent cloud cover and to land versus ocean cover. The clouds, in increments 
from cloud-free to 100% cloud cover, are simulated to be present only on that part 
of the surface covered by vegetation. The violet curve is a simulation with land 
vegetation completely replaced by cloud-free ocean. 

Plankton Essentially all the light absorption that takes place in natural waters is 
attributable to four components of the aquatic ecosystem: the water itself, dissolved 
yellow pigments, the photosynthetic biota (phytoplankton and macrophytes) and 
inanimate particulate matter (tripton) (J. T. O. Kirk, 1983). 

Water absorbs only very weakly in the blue and the green regions of the spectrum, 
but its absorption begins to rise as wavelength increases above 550 nm and is quite 
significant in the red region: a 1-m thick layer of pure water will absorb about 35% 
of incident light of wavelength 680 nm. 

The general tendency is that as phytoplankton concentration increases, reflectance 
decreases in the blue (400-515 nm) and increases in the green (515-600 nm). The 
increased reflectance in the green may be attributed to the fact that phytoplankton, 
being refractive particles, increases scattering at all wavelengths, but in this spectral 
region absorbs only weakly. 

About 1% of the light a photosynthesizing cell absorbs is re-emitted as fluores- 
cence, with a peak at about 685 nm (fig. fTHj) . This fluorescence can show up as 
a distinct peak in the spectral distribution of the upwelling stream or in the curve 
of apparent reflectance against wavelength. Calculations indicate that the increased 
fluorescence associated with an increase in phytoplankton chlorophyll of 1 mg m -3 
in the water would lead to an additional upward radiance of 0.03 W m -2 sr _1 /im" 1 
above the water (Kirk, 1983). 

Plankton also features a red-edge, but it is harder to detect than that for land 
plants. 

In fig. El we show disk averaged spectra of an ocean view of Earth, with an in- 
creasing concentration of plankton. This is an "optimistic" result, since the spectra 
are cloud-free, and we have chosen a particular viewing geometry and a particular 
season, in which plankton can be better detected (NASA-Goddard SeaWiFS web- 
site). The reflectance spectra used as input for our model were taken from (Gower 
et al., 2004) and (Kirk, 1983). 

In fig. El we have increased the concentration of plankton and compared our 
result with the earthshine spectrum. When clouds are taken into account the pres- 
ence plankton is less detectable. In the month of the year when earthshine was 
recorded (June) we do not expect to have too much plankton. Adding the contribu- 
tion of plankton to our simulation does not necessarily improve the fit between the 
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Fi glire 13 ! Fig. at the top: disk- averaged spectra in the optical of a cloud-free Earth seen from 
different views following the diurnal rotation of the planet (ordinates in Rad/Rad top ). We show also 
below examples of the viewing geometries visible during the diurnal rotation (12 am, 6 am, 12 pm 
and 6 pm): as we can see from these pictures the land versus ocean fraction changes considerably 
with time. Fig on the bottom: light curves for a cloud-free, 30% and 60% cloud cover rotating Earth. 
Here we calculate the disk-averaged NDVI of our model using in eq Q) .65-. 68 for the red, .75-. 8 
for the NIR (band I fig. on the left) and .855-. 875 for the NIR (band II fig. on the right). For the 
cloud free case, when the NDVI is less or equal zero, the land fraction covered by vegetation is not 
enough to produce a potentially detectable signal (less than 20%). As we can deduce from the figures 
the band II produces a stronger signal for the clear sky case, but when clouds are included, the NDVI 
threshold is considerably higher (table 2): band II appears less sensitive to the vegetation signal. 
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Fi glire 14: Daily average of the sequence shown in the fig. at the top. In this cloud-free case we 
can still detect the red- edge even if we average over time. 
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Figure 15: Disk-averaged spectra in the optical of an increasingly cloudy Earth. The clouds are 
simulated to be present only on the part of the surface covered by vegetation. From the vantage point 
selected (~ 40 % cloud cover) it is easy to detect the red-edge signal (yellow-green curve, cloud-free). 
As a comparison we have included the same simulation with ocean replacing the land vegetation 
(violet curve). The grey curve at the top shows the situation when clouds are totally covering the 
vegetation. The other plots show intermediate situation (30% and 60% cloud cover). 



data and the model. For high concentrations there is a better agreement with the 
recorded spectrum for wavelengths longer than 0.73 /im, since the red-edge signal is 
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increased, but the fit is worse for shorter wavelengths. For low concentrations the 
contribution of plankton is negligible. 
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Figure 16: Disk-averaged synthetic spectra of a cloud-free, predominantly ocean-view of Earth. 
The concentration of plankton has been increased from (sterile ocean, blue plot) to 300 mg/m 3 
(extreme concentration, red plot). For a present day Earth the average concentration floats between 
0.1 and 10 mg/m 3 , with picks up to 30 mg/m 3 . 




23 



0.05 



0.00 



0.5 



0.6 



0J 



Earthshine 
Model 

4- Plankton 3mg 
+ Plankton 30mg _ 
+ Plankton 300mg 



0.9 



Figure 17: In this figure we show the comparison between the earth-shine spectrum (black plot) and 
our simulation for a sterile ocean (blue plot). We have then progressively increased the concentration 
of plankton in the ocean to test if its presence can improve our fit. 
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Simulation of a TPF detection of an Earth-like planet 

Our synthetic disk-averaged spectra were run through a TPF observation system 
simulator for different spectral resolutions. The TPF book design (TPF book, 1999) 
has been assumed for these calculations and the planet was placed around a G star 
that is lOpc distant. TPF will provide only disk-averaged spectra with possible 
spectral resolving power (X/SX) of ~75 (visible) and ^25 (MIR), depending on the 
final architecture (TPF book), (Beichman and Velusamy, 1999). 




Fi gure 18 ! Simulation of TPF-Interferometer detection of a Earth-like planet orbiting around 
a G star. The top spectrum in each set is at high-resolution, and the middle and lower panels show 
R^lOO (at 10 days integration) and R~20 (at 2 days integration) respectively. 1-a error bars are 
shown in red. We recall that the wavelength range currently proposed and desirable should be 6.-17 
microns for the interferometer, with spectral resolution of 25-50. 
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Fi glire 19l Simulation of TPF-Coronograph detection of a Earth-like planet orbiting around a G 
star. These plots have been produced by processing the synthetic spectra with an observational system 
simulator of the TPF coronograph, at spectral resolutions R= 70 with corresponding integration time 
of 6h. The error bars in the spectra represent rms noise (1 sigma) in the TPF measurement in each 
spectral channel for the integration. 
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5 Discussion 



Ford et al. (2002), as we already mentioned in the introduction, developed a model 
to evaluate the photometric variability in some specific bands in the optical (450, 
550, 650, and 750 nm) due to daily rotation and seasons. In their paper, they 
predict a daily variation as high as 150% for their cloud- free Earth model (750 nm 
band). According to our results for the clear sky case (see fig. [T3|) . the maximum 
daily variation is about 50% (750 nm band). The seasonal variation though, can be 
considerably higher due to the changes in ice cover, which can dramatically increase 
the surface albedo (fig. IHIand UHjl . Clouds tend to increase the overall brightness and 
variability, reducing the fractional variation in the reflected light curve compared to 
the cloud free case (fig. ITTTI and fT5|) . These conclusions although, in the case 
of an extrasolar terrestrial planet, should be taken with some caveats, since they 
depend on the cloud pattern and rate of cloud formation on that planet. We confirm 
their hypothesis that the diurnal rotation variation of an Earth-like planet is lower 
in the IR than in the optical (fig. fT2|) . since the temperature gradients are not so 
dramatic compared with the surface albedo variations. This is due to the presence of 
a relatively thick atmosphere and of an ocean that act as buffer to the temperature 
excursions that otherwise would occur (see for comparison the Mars case (Tinetti et 
al, 2004)). 

Previous workers have utilized Earthshine data to look for a vegetation signature 
in a disk-averaged spectrum (Woolf et al., 2002; Arnold, et al., 2002). In the case 
of Arnold et al., we cannot directly compare our results with theirs, since we do not 
know the details of these observations (spectra recorded, viewing geometries, time, 
illumination, etc.). We agree, though, at least qualitatively, with their conclusions. 

The Earthshine data show that the atmosphere on Earth is clear enough to allow 
some of the vegetation signal to show through, but it is weakened. Our model and 
disk- and diurnally-averaged NDVI show that the vegetation signal has the potential 
to be quite strong but also completely hidden (fig. Q31 El an d fl5|) . The problem 
of clouds obscuring vegetation on land is compounded by the fact that vegetation 
promote conductance of moisture from the soil to the atmosphere, and vegetation 
roughness promotes moisture convergence and hence cloud cover at the regional and 
continental scales (Dickinson, 1989; Friend and Kiang, accepted). Cloud cover will 
tend to favor vegetated areas (and vice versa). The role of transpiration processes 
in the hydrological cycle on an extrasolar planet, should therefore be considered 
in narrowing the range of plausible detectable biosignatures. Choice of bands for 
calculation of the red edge detectability must clearly take into account the obscuring 
of surface signatures not only by clouds but also by other atmospheric gases. Oxygen 
is a potential biosignature produced by photosynthesis and water is a prerequisite for 
known life; therefore, both will be present if a vegetation red edge is present, while 
both obscure parts of the red edge. The choice of red and NIR bands for our red 
edge index (disk- average NDVI) followed the choices of .65-. 68 for the red and .75-. 8 
(band I) or .855-. 875 (band II) for the NIR since the O2 and H2O absorbance lines 
all fall along the red edge feature. NIR band I is subject to physiological variability, 
whereas the band II wavelength range is fairly stable, even across mosses, lichens, and 
cyanobacteria. Unfortunately NDVI band II loses its sensitivity when we consider 
the impact of level of cloud cover on the strength of the vegetation signature (fig. 

With spectral resolutions of ~ 75 in the visible for TPF-C, it will be important 
that the ranges of the sensor bands be designed both to detect atmospheric gases, 
but also to exclude them for surface sensing. Given constraints of sensor design 
and damping effects of clouds and diurnal averaging, we will need to determine how 
strong an NDVI signal is an actual indicator of photosynthesis. 
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Therefore, unless we uncover a potentially abiotic source of a red edge, we may be 
able to say that any statistically significant departure of the NDVI from zero implies 
the existence of Earth-like land-based photosynthetic activity. Given potential TPF- 
C sensor bands at 650-680 nm and 750-800 nm, the red-edge detectability for the 
Earth requires at least 20% average diurnal land vegetation cover unobscured by 
clouds to achieve an index that is greater than zero. So, this means a planet with a 
high percent of water over the surface could still reveal land surface photosynthetic 
activity. When clouds are included, the thresholds for discriminating unambiguously 
the red-edge become 0.03 (band I) and 0.08 (band II). These numbers although might 
be significantly smaller and even negative depending on phase, viewing geometries 
and cloud distribution. 

Moving from land to marine vegetation we have dedicated a paragraph on the 
detectability of plankton (fig. [TH| IT7|> . In fig. QSlwe show disk averaged spectra of 
an ocean view of Earth, with an increasing concentration of Plankton. As we have 
already pointed out, this is an "optimistic" result, since we have supposed to have a 
clear sky, and we have chosen a particular viewing geometry and a particular season, 
in which Plankton can be better detected (NASA-Goddard SeaWiFS website). When 
clouds are taken into account, a more realistic case, (fig. IT7|> plankton is very hard 
to detect, even for high concentration. 

The concentration of 300 mg of Chlorophyll is not realistic for our present day 
planet, but it is useful to understand the trend of how plankton can possibly modify 
the final spectra. In the case of an early Earth or an extrasolar terrestrial planet this 
hypothesis might be plausible. It is out of the interests of this paper to discuss further 
the likelihood of the occurrence of photosynthesis on Earth-like planets outside our 
solar system on which life has evolved. Based on knowledge of photosynthesis on 
Earth and of stellar evolution Wolstencroft and Raven in (Wolstencraft and Raven, 
2002) conclude that it is likely that photosynthesis would have evolved on Earth- 
like planets in response to the same evolutionary factors as have been involved on 
Earth. The chemical intermediates and catalysts would probably be different, but the 
substrates and end-products would have been the same. Modeling of photosynthesis 
on Earth-like planets orbiting stars of different spectral types shows that cooler 
stars, with maximum radiation output at longer wavelengths, may require more 
than the two light reactions used in oxygen-evolving photosynthesis on Earth; such 
photosynthesis would be limited by the attenuation of radiation by water. Detection 
of photosynthesis will be based on spectroscopy using most probably photosynthetic 
pigments which could have very different absorption properties from those on Earth. 
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6 Conclusions 



• Our approach is feasible. The comparison with the experiment confirms: the 
model works (see fig. |HJ)! 

• We can distinguish among different surface types looking at the Earth's spectra 
in the visible (0.7-0.8 fim band, red-edge), there are almost no differences in 
the IR (fig. 

• The contribution of clouds is dramatic. We have to include them in order to 
have a realistic model (see fig. El El EJ 

• Our disk-averaged spectra are sensitive to phases above all in the optical and 
in particular when clouds are present (fig. fTTl). 

• We have generated disk- aver aged spectra from different vantage points. In the 
optical the main differences are due to a changed surface composition. In the 
IR on the contrary the fluctuations are due to the temperature gradients. 

• The lack of high spectral resolution smears out the spectral features. However 
some of the strongest features are detectable with TPF (fig. ITH| fTHj) . 

• Time dependent variations in the disk-averaged spectra, or "light-curve" can 
provide additional information about spatial variations (fig. E1E1EJ)- 

• The red-edge biosignature could be potentially detectable (depending on the 
instrumentation capabilities) when we consider a cloud-free Earth. If clouds 
are randomly distributed and in high %, it is very hard to detect this signal 

(fig. ESI). 

• Plankton could be potentially discriminated in disk-averaged spectra, depend- 
ing on the concentration, season of the year, viewing geometry and above all 
cloud cover (fig. El and EI)- 
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